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A brief review is given of recent results devoted to the effects of large-scale anisotropy on the 
inertial-range statistics of the passive scalar quantity #(£, x), advected by the synthetic turbulent 
velocity field with the covariance oc S(t — t')l x — x'| e . Inertial-range anomalous scaling behavior 
is established, and explicit asymptotic expressions for the structure functions S n (r) = ([0(t, x + 
r) — #(t,x)] n ) are obtained; they are represented by superpositions of power laws with universal 
(independent of the anisotropy parameters) anomalous exponents, calculated to the first order in 
e in any space dimension. The exponents are associated with tensor composite operators built of 
the scalar gradients, and exhibit a kind of hierarchy related to the degree of anisotropy: the less 
is the rank, the less is the dimension and, consequently, the more important is the contribution to 
the inertial-range behavior. The leading terms of the even (odd) structure functions are given by 
the scalar (vector) operators. The small-scale anisotropy reveals itself in odd correlation functions: 
for the incompressible velocity field, Sz/S 2 decreases going down towards to the depth of the 
inertial range, while the higher-order odd ratios increase; if the compressibility is strong enough, the 
skewness factor also becomes increasing. 
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> 
ly-s ' The investigation of intermittency and anomalous scaling in fully developed turbulence remains essentially an open 

theoretical problem. Both the natural and numerical experiments suggest that the deviation from the predictions of 

the classical Kolmogorov-Obukhov theory! R0 is even more strongly pronounced for a passively advected scalar field 

than for the velocity field itself; see, e.g., [|3||| and literature cited therein. At the same time, the problem of passive 

advection appears to be easier tractable theoretically: even simplified models describing the advection by "synthetic" 

velocity fields with given Gaussian statistics reproduce many of the anomalous features of genuine turbulent heat 

or mass transport observed in experiments. Hence the term "toy models." Therefore, the problem of passive scalar 

advection, being of practical importance in itself, may also be viewed as a starting point in studying anomalous scaling 

-H in the turbulence on the whole. 

The most remarkable progress has been achieved for the so-called rapid-change model, which dates back to Obukhov 
^- , and Batchelor and was explicitly introduced by Kraichnan |5|: for the first time, the anomalous exponents have been 
calculated on the basis of a microscopic model and within regular perturbation expansions; see, e.g., [6-10] and 
references therein. 

Within the "zero- mode approach," developed in [6-9] , nontrivial anomalous exponents are related to the zero modes 
(unforced solutions) of the closed exact equations satisfied by the equal-time correlations. Within the approach based 
on the field theoretic renormalization group (RG) and operator product expansion (OPE), the anomalous scaling 
emerges as a consequence of the existence in the model of composite operators with negative critical dimensions, 
which determine the anomalous exponents [10-13]. 

An important question recently addressed is the effects of large-scale anisotropy on the inertial-range statistics of 
passively advected fields [8,9,12-18] and the velocity field itself 

According to the classical Kolmogorov-Obukhov theory 
forcing (boundary conditions, geometry of an obstacle etc) dies out when the energy is transferred down to smaller 
scales owing to the cascade mechanism. A number of recent works confirms this picture for the even correlation 
functions, thus giving some quantitative support to the aforementioned hypothesis on the restored local isotropy 
restoration of the inertial-range turbulence for the velocity and passive fields [12-19]. More precisely, the exponents 
describing the inertial-range scaling exhibit universality and hierarchy related to the degree of anisotropy, and the 
leading contribution to an even function is given by the exponent from the isotropic shell [12,15-19]. Nevertheless, 
the anisotropy survives in the inertial range and reveals itself in odd correlation functions, in disagreement with what 
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the anisotropy introduced at large scales by the 



was expected on the basis of the cascade ideas. The so-called skewness factor decreases down the scales much slower 
than expected |§j4|g,g|, while the higher-order odd dimensionless ratios (hyperskewness etc) increase, thus signalling 
of persistent small-scale anisotropy |I^,|lJ,|lG| . The effect seems rather universal, being observed for the scalar [|12J and 
vector jl(| fields, advected by the Gaussian rapid-change velocity, and for the scalar advected by the two-dimensional 
Navier-Stokes (NS) velocity field @. 

Below we briefly review the most recent results obtained on the example of the rapid-change model using the RG 
and OPE techniques; more detailed discussion can be found in the original papers [10-13]. The detailed description 
of the RG approach to the models of fully developed turbulence is given in the book [£(J and review paper |2l]] . For 
simplicity, we restrict ourselves with the rapid-change models; the case of a finite-correlated velocity field is considered 

n& 

The advection of a passive scalar field in the presence of an imposed linear gradient is described by the equation 

d t 6 + (vd)0 = v d 2 9 - (hv). (1) 

Here 9(x) = 0(£, x) is the random (fluctuation) part of the total scalar field Q(x) = 6(x) + (hx), h is a constant 
vector that determines distinguished direction, vq is the molecular diffusivity coefficient, dt = d/dt, di = d/dxi, and 
d 2 = didi is the Laplace operator. 

The velocity field v(x) = {vi(x)} obeys a Gaussian distribution with zero mean and correlator 

(v i (x)v j (x')) = S(t-t')K ij (r) (2) 

with 

^ (r) = Z3o /(2^W eXp[ik ' r] ' r " X ~ X '' (3) 

where Py(k) = <% — kikjk 2 is the transverse projector, k is the momentum, k = |k|, Dq > is an amplitude factor, 
and < e < 2 is a free parameter. The IR regularization is provided by the cut-off in the integral (0) from below 
at k a m, where m = l/L is the reciprocal of the integral turbulence scale; the precise form of the cut-off is not 
essential. For < e < 2, the difference 



has a finite limit for m — ► 0: 



with 



S l0 (r) = K l0 (S)) - K %3 {v) (4) 



S l0 {r)=Dr £ (d + e-l)S i:j -e-^- , (5) 



D _ -D T(-e/2) 



(47r) d / 2 2^(d + e)r(d/2 + e/2) ' 

where T(- • •) is the Euler gamma function (note that D > 0). It follows from Eq. (0) that e can be viewed as a kind 
of Holder exponent, which measures the roughness of the velocity field. In the RG approach, the exponent e plays 
the same role as the parameter e = 4 — d does in the RG theory of critical phenomena (P3,p3[ . The relations 

g a = D /u = A E (6) 

define the coupling constant go (i.e., the expansion parameter in the ordinary perturbation theory) and the charac- 
teristic ultraviolet (UV) momentum scale A. 

The term (hv) in Eq. (Q) acts as a driving force; it maintains the steady state and, simultaneously, introduces 
the large-scale anisotropy. In the original (isotropic) formulation of the model this term is replaced with an artificial 
Gaussial random force f(x) with the covariance 

(f(x)f(x')) = 6(t-t')C(mr), r=|x-x'|, (7) 

where C{mr) is a function finite as mr = 0. 

The issue of interest is, in particular, the behavior of the equal-time structure functions 

S„(r) = <[0(i,x)-0(t,x')n (8) 



in the inertial range, specified by the inequalities A >> 1/r >> to. The dimensionality considerations give 

S n (r) = v- n/2 r n R n (Ar 7 mr,z), (9) 

where z is the cosine between the directions of r and h and R n are some functions of dimcnsionlcss parameters. 
In principle, they can be calculated within the ordinary perturbation theory (i.e., as series in go), but this is not 
useful for studying inertial-range behavior: the coefficients are singular in the limits Ar — > oo and/or mr — * 0, which 
compensate the smallness of go (assumed in ordinary perturbation theory), and in order to find correct IR behavior 
we have to sum the entire series. The desired summation can be accomplished using the RG and OPE; see Refs. 

us 

The RG analysis consists of two main stages. On the first stage, the multiplicative renormalizability of the model is 
demonstrated and the differential RG equations for its correlation functions are obtained. The asymptotic behavior 
of the functions like (0) for Ar > > f and any fixed mr is given by IR stable fixed points of the RG equations and has 
the form 

S n (r) = D~ n/2 r n ^W R n (mr, z), Ar » 1, (10) 

with certain, as yet unknown, scaling functions R n (mr, z) . 

On the second stage, the small mr behavior of the functions R n is studied within the general representation (nfl) 
using the OPE. It shows that, in the limit mr — > 0, the functions R n have the asymptotic forms 

R n (mr 7 z) — 2^CF{mr, z) (mr) F , (11) 

F 

where Cf are coefficients regular in mr. In general, the summation is implied over all possible renormalized composite 
operators F allowed by the symmetry; Ap being their critical dimensions. 

The peculiarity of the models describing turbulence is the existence of the so-called "dangerous" composite operators 
with negative critical dimensions |ljj [12],^0 21 1 . Their contributions into the OPE give rise to a singular behavior of 



the scaling functions for mr — > 0, and the leading term is given by the operator with minimal A^ . In the isotropic 



case only scalar operators have nonvanishing mean values and contribute to the right hand side of Eq. (11). The odd 
functions S2/C+1 vanish, while the leading contributions to S 2 k are determined by scalar gradients Fk — (di8di6) h and 
have the form 

S 2k (r) ~ V > fe(2 ~ E) (™-) A, % (12) 

where the critical dimensions A^ of the operators Fk are given by 

A fe = -2fc(fc - l)e/(d + 2) + 0(e 2 ) = -2k(k - l)e/d + 0(1/ d 2 ). (13) 



The expression ( |l2| ) agrees with the results obtained earlier in |]|7|] using the zero-mode techniques; the 0(e 2 ) contri- 
bution to Afc is reported in JlQ]. Note that the operators Fi with I > k (whose contributions would be more important) 
do not appear in the right hand side of the OPE for S 2 k ', this is a consequence of the linearity of the original equation 
@m6(x). 

In the anisotropic case, an important role is played by the irreducible tensor operators of the form 

F[n, P ] = d n e--- d lp e (diOdie) 1 + ■ • • , (w) 

where p is the number of the free vector indices and n = p + 21 is the total number of the fields 9 entering into 
the operator; the vector indices of the symbol F[n,p] are omitted. The dots stand for the subtractions (needed 
for p > 2) which ensure that the resulting expression is traceless with respect to any contraction, for example, 
F[2, 2} = drfdjO - SijFft, 0]/d. 

For h ^ 0, the operators (|14| ) acquire nonvanishing mean values, built of the vectors h and delta symbols, and give 
nontrivial contributions to the OPE for S n ; in particular, this gives rise to nonvanishing odd functions. 

The critical dimensions associated with the operators ( |l4| ) have the form [Q (see also Refs. |l^ , p4| ): 

A[nri = ;h ,(-l)-(^l)(.-rtg + n + y -2) t + 0(t3) (15) 

The straightforward analysis of the explicit one-loop expression ([15]) shows that for any fixed n, the dimension 
A[ra,p] decreases monotonically with p and reaches its minimum for the minimal possible value of p — p n , i.e., p = 
if n is even and p = 1 if n is odd: 



A[n,p] > A[n,p'] if p > p'. (16) 

Furthermore, this minimal value A[n,p n ] is negative and decreases monotonically as n increases: 

> A[2fc, 0] > A[2k + 1, 1] > A[2fc + 2, 0]. (17) 

Finally, we note that for any fixed p, the dimension (|15|) decreases monotonically as n increases: 

A[n,p] > A[n',p] if n<ri. (18) 

The mean value of the operator with the dimension A[fc,p] is a traceless (for p > 2) tensor built of the vector h 
and Kronecker delta symbols; the contraction with the vector indices of the corresponding coefficient Cf in Eq. (ThD) 
gives rise to p-th order Legendre polynomial P p (z) 7 so that the expansion ([ll]) is analogous to the decomposition in 
irreducible representations of the rotation group, used in Refs. [15-19]. 

The inequalities (fig), ( |17| ) show that the contributions of the tensor operators (Q) into the asymptotic expression 
( |ll| ) exhibit a kind of hierarchy: the less is the rank, the more important is the contribution. 

The leading term of the expression ([ll]) for the even (odd) function S n is determined by the scalar (vector) composite 
operator consisting of n factors d0 and has the form 

S n (r) oc (hr) n K) A I" A ' . (19) 

Let us now explore the consequences of our results for the issue of the local isotropy restoration in this problem. 
As a measure of small-scale anisotropy one usually uses the dimcnsionless ratios of the structure functions: 

R n (v)=S n (r)/S™ /2 (v). (20) 

From Eq. (|l^) it then follows that in inertial range we have: 

R 2k+1 OC ( mr )A[2fc+M]-(2fc+l)A[2,0]/2 ; pi) 

i? 2fe oc(mr) A [ 2fc ^- fcA ^°]. (22) 

The dependence on the Peclet number, Pe = (A/m) £ can be estimated by replacing r with 1/A; see Ref. |§. Using 
explicit O(e) expressions for A[n,p] we then obtain: 

i? 2fe+1 ocP e -( d+2 - 4fe2 )/[ 2 ( d + 2 )], (23) 

R 2k ozPe 2k{k - l ^ {d+2 \ (24) 

Since the leading terms of the even functions (n9h are determined by the exponents of the isotropic shell (i.e., those 
related to the scalar composite operators), the inertial-range behavior of the even ratios J22]), J24| ) is the same as in 
the isotropic model. This gives a quantitative support to the universality of leading anomalous exponents for the even 
functions with respect to different classes of forcing. On the other hand, the odd quantities (pi|), (23|) appear sensitive 
to the anisotropy: R3 in (J23T) slowly decreases for Pe — > 00, while ratios i?2fe+i with k > 2 increase with Pe, thus 
signalling of the persistent small-scale anisotropy. 

The picture outlined above seems rather general. Indeed, it is compatible with that recently established for the 
NS turbulence |19| , the scalar field passively advected by a NS velocity in the two-dimensional inverse energy cascade 
ipTf , and for the vector (magnetic) field advected by the velocity of the type (||), see Refs. |15Hl7|. 

To conclude with, let us briefly discuss the effects of compressibility onto the small-scale anisotropy persistence. 
The transversal projector in the correlator (|3|) is then replaced by the combination 

P 4i (k) + aQ«00. 

where Qy(k) = kikj/k 2 is the longitudinal projector and a > is a free parameter. The model (nT) then corresponds 
to the advection of a tracer (entropy, temperature, concentration of a pollutant), while the nonlinearity Q'd{v6) 
describes advection of a density. 

Equation (|l5|) for a tracer becomes |13| 



1 _ 2n(n - l)(l - g) - ( n - p)(n + p + d - 2)(d + 1 + a) 
A[n ' Pl 2(d+2)(d-l + a) £ ' (25j 

the 0(e 2 ) correction is calculated in fl25|| for all n, p, d and a. 

Although the hierarchy relations fll6|)-(|l8|) remain valid for all values of a > (in particular, dA[n : p]/dp > 0), the 
corrections become closer to leading terms as a increases: 

d 2 A[n, P ]/dpda<0. 



Furthermore, the expressions (23) take on the form 



R 2k+1 oc (jnrf l(d-l+ a )(d+2-4k 2 )-8a k 2 ]/2(d+2)(d-l+ a ) _ /ggs 

It follows from ( p6| ) that -R3 diverges for mr — » provided a is large enough [namely, a > (d— l)(d+2)/(10— d)+0(e)], 
while the divergence of the higher ratios i?2fc+i enhances as a increases. 

This means that compressibility enhances the penetration of the large-scale anisotropy towards the depth of the 
inertial range. This fact also seems universal, being observed in the model of the passively advected magnetic field 

II- 

Some remarks are now in order. 

The critical dimensions of all composite operators ( Jl 5|) are independent of the forcing: they remain unchanged, 
when the isotropic stirring force rt7h is replaced in Eq. (pi) with the anisotropic term (hv). The difference is that 
for the former case, only scalar operators contribute to the representations like (Jll|) , while for the latter case the 
irreducible tensor operators acquire nonzero mean values and also contribute to Eq. (|ll|). This is easily understood 
in the language of the zero-mode approach: the noise / and the term (hv) do not affect the differential operators in 
the equations satisfied by the equal-time correlations functions; the zero modes (homogeneous solutions) coincide in 
the two cases, but in the latter case the modes with nontrivial angular dependence should be taken into account. In 
the language of the RG (which is also applicable to the case of a finite-correlated velocity, see Refs. |1J,[r|) this is 
explained as follows: the stirring force or field h do not enter into the diagrams that determine the renormalization 
of the operators ( |l5|) . However, the anomalous exponents of the passive scalar become nonuniversal and acquire the 
dependence on the anisotropy parameters if the velocity field is taken to be strongly anisotropic also at small scales 

The picture outlined above and in Refs. |l5],[l6| for passively advected fields (a superposition of power laws with 
universal exponents and nonuniversal amplitudes) seems rather general, being compatible with that established re- 
cently in the field of NS turbulence, on the basis of numerical simulations of channel flows and experiments in the 
atmospheric surface layer, see Refs. |19| and references therein. In those papers, the velocity structure functions 
were decomposed into the irreducible representations of the rotation group. It was shown that in each sector of the 
decomposition, scaling behaviour can be found with apparently universal exponents. The amplitudes of the various 
contributions are nonuniversal, through the dependence on the position in the flow, the local degree of anisotropy and 
inhomogeneity, and so on Jl9[ . 

This is rather surprising because the equations for the correlation functions in such cases are neither closed nor 
isotropic and homogeneous. Although the hierarchy similar to Eq . (fij^ ) is demonstrated by the critical dimensions 
of certain tensor operators in the stirred NS turbulence, see Ref. [j28| and Sec. 2.3 of |2y], the relationship between 
them and the anomalous exponents is not obvious there. One can thus speculate that the anomalous scaling for the 
genuine turbulence can also appear a linear phenomenon. It is worth recalling here that the so-called "additive fusion 



rules," hypothesized for the NS turbulence in a number of papers, Refs. |29 , and characteristic of the models with 
multifractal behavior (see Ref. J30||), arise naturally in the context of the rapid-change models owing to their linearity. 
The existing results for the Burgers turbulence can also be interpreted naturally as a consequence of similar fusion 
rules, where only finite number of dangerous operators contributes to each structure function, see Ref. |3|] . 

Of course, one should not insist too much on this bold assumption. 

To conclude, let us compare briefly the situation for the passively advected fields with the case of weak acoustic 
turbulence, where the spectra can be obtained as solutions of the linear kinetic equations (see Refs. |3^,|33|]). For weakly 
dispersive waves (e.g., with the dispersion law u)(k) ex k 1+ with S << 1), the anisotropy introduced by the large-scale 
forcing enhances going down towards to the depth of the inertial range [ p2| . The hierarchy of the exponents related 
to the Legendre decomposition is opposite to that established below and in Refs. |l2|Jl^ , pj| : anisotropic corrections 
decrease slower for larger p's [[52|. On the contrary, for the nondispersive waves (5 = 0) the hierarchy of the exponents 
is similar to that in our case, the anisotropic corrections decay faster and faster with p and the spectrum tends to 
become isotropic at small scales J3JJ. To the best of our knowledge, no information is available for the higher-order 
correlation functions for such models. One can thus conclude that turbulent systems can exhibit essentially different 
types of behavior with respect to the small-scale isotropy restoration. 
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